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Abstract. The Lunar Laser Ranging (LLR) experiment provides precise observations 
of the lunar orbit that contribute to a wide range of science investigations. In 
On ' particular, time series of highly accurate measurements of the distance between the 

Earth and Moon provide unique information that determine whether, in accordance 
O ' with the Equivalence Principle (EP), both of these celestial bodies are falling towards 

^ , the Sun at the same rate, despite their different masses, compositions, and gravitational 

Jh ' self-energies. Analyses of precise laser ranges to the Moon continue to provide 

increasingly stringent limits on any violation of the EP. Current LLR solutions give 
(—0.8 ± 1.3) X 10~^^ for any possible inequality in the ratios of the gravitational and 



inertial masses for the Earth and Moon, {mQ/'mi)-E, — {mQ/'mi)M- Such an accurate 
C^ . result allows other tests of gravitational theories. Focusing on the tests of the EP, we 

ly^ I discuss the existing data and data analysis techniques. The robustness of the LLR 

solutions is demonstrated with several different approaches to solutions. Additional 
high accuracy ranges and improvements in the LLR data analysis model will further 
advance the research of relativistic gravity in the solar system, and will continue to 



o 

p^ i provide highly accurate tests of the Equivalence Principle 
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1. Introduction 

The Equivalence Principle (EP) lies at the foundation of Einstein's general theory of 
relativity; testing this fundamental assumption with the highest possible sensitivity is 
clearly important, particularly since it may be expected that the EP will not hold in 
quantum theories of gravity [25]. In general relativity, test objects in free-fall follow 
geodesies of space-time, and what we perceive as the force of gravity is instead a result 
of our being unable to follow those geodesies of space-time, because the mechanical 
resistance of matter prevents us from doing so. One can test both the validity of the 
EP and of the field equations that determine the geometric structure created by a 
mass distribution. There are two different "flavors" of the Principle, the weak and the 
strong forms of the EP that are currently tested in various experiments performed with 
laboratory test masses and with bodies of astronomical sizes. 

In its weak form (the WEP), the Principle states that the gravitational properties 
of primarily strong and electro-weak interactions obey the EP. In this case the relevant 
test-body differences are their fractional nuclear-binding differences, their neutron-to- 
proton ratios, their atomic charges, etc. General relativity and other metric theories of 
gravity assume that the WEP is exact. However, many gravitational theories predict 
observable violations of the EP at various fractional levels ranging from 10~^^ down to 
10~^^ [6]. For instance, extensions of the Standard Model (for discussion, see ^M 125]) 
that contain new macroscopic-range quantum fields predict quantum exchange forces 
that generically violate the WEP because they couple to generalized "charges" rather 
than to mass/energy alone, as in general relativity. Therefore, even a confirmation that 
the WEP is not violated at some level will be exceptionally valuable, placing useful 
constraints on the range of possibilities in the development of a unified physical theory. 

In a laboratory, precise tests of the EP can be made by comparing the free fall 
accelerations, ai and 02, of different test bodies. When the bodies are at the same 
distance from the source of the gravity, the expression for the EP takes the form 

Aa 2(ai — 02) 



.mj. 


1 


-tug' 
.mj. 


= A 
2 


'mo' 
.mi . 



(1) 

a (01-1-02) Lm/Ji Vmj\2 Lm/J 

where niQ and m/ are the gravitational and inertial masses of each body. The sensitivity 
of the EP test is determined by the precision of the differential acceleration measurement 
divided by the degree to which the test bodies differ (e.g. composition). 

Currently, the most accurate results in tests of the composition-independence of 
acceleration rates of various masses toward the Earth were reported by ground-based 
laboratories p]. A recent experiment measured the fractional differential acceleration 
between Be and Ti test bodies at the level of Aa/a = (0.3 ± 1.8) x 10"^^ |21]. The 
accuracy of these experiments is high enough to confirm that the strong, weak, and 
electromagnetic interactions each contribute equally to the passive gravitational and 
inertial masses of the laboratory bodies. A review of the recent laboratory tests of 
gravity can be found in [1] . Significant improvements in the tests of the EP are expected 
from dedicated space-based experiments [25] . 
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In its strong form (the SEP) the EP is extended to cover the gravitational properties 
resulting from gravitational energy itself, thus involving an assumption about the non- 
linear property of gravitation [30]. In the SEP case, the relevant test body differences 
are the fractional contributions to their masses by gravitational self-energy. Although 
general relativity assumes that the SEP is exact, many modern theories of gravity 
typically violate the SEP by including new fields of matter, notably scalar fields [3], H]. 
The ratio of the gravitational-to-inertial masses for a body can be parameterized as 

where m is the mass of a body, U is the body's gravitational self-energy (t/ < 0), mc^ 
is its total mass-energy, and r^ is a dimensionless constant for SEP violation (77 = in 
general relativity, see discussion in [3D]). Because of the extreme weakness of gravity, a 
test of the SEP requires bodies of astronomical sizes. 

At present, the Earth- Moon-Sun system provides the best solar system arena for 
testing the SEP. Time series of the highly accurate measurements of the distance 
between the Earth and Moon provide unique information used to determine whether, in 
accordance with the EP, both of these celestial bodies are faUing towards the Sun at the 
same rate, despite their different masses, compositions, and gravitational self-energies. 

This paper focuses on the tests of the EP with LLR. To that extent. Section [2] 
discusses the LLR history, experimental technique, and the current state of the effort. 
Section [3] discussed LLR range data and distribution. Section H] focuses on the precision 
tests of the EP with LLR data and discusses the results obtained. Section |5] discusses 
future prospects. In Section [6] we conclude with a summary and outlook. 

2. LLR Technique 

Over 40 years since their initiation, analyses of precision laser ranges to the Moon 
continue to provide increasingly stringent limits on any violation of the EP. 

Each LLR measurement is the round-trip travel time of a laser pulse between an 
observatory on the Earth and one of the five corner-cube retro-reflector (CCR) arrays 
on the moon. To range the moon, the observatories fire a short laser pulse toward 
the target array. The lasers currently used for ranging operate at 20 Hz, with a pulse 
width of about 200 psec; each pulse contains ~ 10^® photons. Under favorable observing 
conditions a single reflected photon is detected every few seconds for most LLR stations 
and in less than one second for Apache Point [lOj. Such a low return rate is due to the 
huge attenuation encountered during the round-trip of the pulse. The outgoing narrow 
laser beam must be accurately pointed at the target since the beam's angular spread is 
typically a few arcsec; the spot size on the moon is a few km across. The amount of 
energy falling on the CCR depends inversely on that spot area. 

The returning pulse illuminates an area around the observatory that is a few tens 
of kilometers in diameter (~10 km for the 532 nm green light and Apollo CCRs). The 
observatory has a sensitive detector which records single photon arrivals. The power 
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received by the telescope depends directly on the telescope's collecting area and inversely 
on the returning spot area. Velocity-caused aberration of the returning beam is (~1 
arcsec. 



2.1. Equivalence Principle and the Earth- Moon system 

The Jet Propulsion Laboratory's (JPL) program which integrates the orbits of the Moon 
and planets considers accelerations due to Newtonian, geophysical and post-Newtonian 
effects. The dynamics of the three-body Sun-Earth-Moon system in the solar system 
barycentric inertial frame provides the main LLR sensitivity for a possible violation of 
the equivalence principle. In this frame, the quasi-Newtonian acceleration of the Moon 
with respect to the Earth, a = slm — sle, is calculated to be: 



a = — /i 
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where /i* = lJi-E{jnG/mj)y[ + fiM{fnG/''T^i)E and /i^ = Gmk. The first term on the right- 
hand side of the equation above, is the acceleration between the Earth and Moon with 
the remaining pair being the tidal acceleration expression due to the solar gravity. The 
above acceleration is useful for either the weak or strong forms of the EP. 
Rearranging the above equation emphasizes the EP terms: 
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The presence of the EP parameters for gravitational-to-inertial mass ratios in /i* modifies 
Kepler's third law to n'^a^ = fi* for the relation between semimajor axis a and mean 
motion n in the elliptical orbit approximation. This term is notable, but in the LLR 
solutions /iE + /iM is a solution parameter with uncertainty, so this term does not provide 
a sensitive test of the EP, though its effect is implicit in the LLR solutions. The second 
term on the right-hand side with the differential acceleration toward the Sun is the 
Newtonian tidal acceleration. The third term involving the mG/mi ratios for two bodies 
gives the main sensitivity of the LLR test of the EP. Since the distance to the Sun is 
~390 times the distance between the Earth and Moon, the last term is approximately 
the difference of any EP violation of the two bodies times the Sun's acceleration of the 
Earth- Moon center of mass. 

Treating the EP related tidal term as a perturbation [TB] found a polarization of 
the Moon's orbit in the direction of the Sun with a radial perturbation 

Ar = S[{mG/mi)E - {mG/mj)M] cos D, (5) 

where S" is a scaling factor of about —2.9 x 10^^ mm (see [161 E])- 

The Earth and Moon are large enough to have significant gravitational self-energies 
and a lunar test of the EP was proposed by Nordtvedt [15] . Both bodies have differences 
in their compositions and self-energies and the Sun provides the external gravitational 
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acceleration. For the SEP effect on the Moon's position with respect to the Earth it is 
the difference of the two self-energy values which is of interest [26] : 

^ ^ ^ ^ ^ 4.45 X 10-1°. (6) 



.mc^/E Kmc^ J M 
For the SEP, combining Eqs. ([2]) and ([5]) yields 

Ar = St] [(u/mc^)^ - (u/mc^)^] cosD, (7) 

Ar = CoT]cosD. (8) 

Applying the difference in numerical values for self-energy for the Earth and Moon 
gives a value of Cq of about 13 m (see pTf [5]). In general relativity 77 = 0. A unit value 
for 7] would produce a displacement of the lunar orbit about the Earth, causing a 13 m 
monthly range modulation. (See Sec. [Hand Table [3] for a comparison of the theoretical 
values of S and Co with numerical results.) 

In essence, LLR tests of the EP compare the free-fall accelerations of the Earth and 
Moon toward the Sun. Lunar laser-ranging measures the time-of-fiight of a laser pulse 
fired from an observatory on the Earth, bounced off of a retro-refiector on the Moon, 
and returned to the observatory (see [3 |30]). If the EP is violated, the lunar orbit will 
be displaced along the Earth-Sun line, producing a range signature having a 29.53 day 
synodic period (different from the lunar orbit period of 27 days). For a review of history 
and present capabilities of LLR to test the EP, see [30] . 



2.2. LLR Normal Points, Model and Science Outcome 

A normal point results from a statistical combination of the observed transit times of 
several individual photons arriving at the observing detector within a relatively short 
time span, typically minutes to tens of minutes [SOj [11] . The resulting "range" normal 
point is the round trip flight time for a particular firing time. In addition to range 
and time, supporting information includes uncertainty, signal to noise ratio, number of 
photons, atmospheric pressure, temperature, and humidity. 

The existing model formulation at JPL and its computational realization in 
computer code, is the product of many years of effort [H [26]. The successful analysis 
of LLR data requires attention to geophysical and rotational effects for the Earth and 
the Moon in addition to orbital effects. The main features of the current model are 
summarized in [26l [28l [22] . For a general review of LLR see [7j. Refs. [27l [28l [29l [30] 
have the most recent results for gravitational physics. 

3. Range Data and Distribution 

The solutions presented here use 17,580 laser ranges from March 1970 to July 2011. 
The ranging stations include the McDonald Observatory in Texas (37.7% of the total 
observations), the Observatoire de la Cote d'Azur (OCA) in France (52.9%), the 
Haleakala Observatory in Hawaii (3.9%), and the APOLLO facility at the Apache 
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Figure 1. (a) Annual rms residuals of LLR data from 1970 to 2011. (b) A histogram 
of the distribution of LLR observations vs. elongation angle D for all of the data. 



Point Observatory in New Mexico (5.4%). The OCA station is described by [20]. The 
APOLLO effort is described in the article by pjy in this issue. The Matera station in 
Italy has demonstrated lunar capability with 0.1% of the total observations, but does 
not range the Moon regularly. There are five retro-reflectors on the Moon at the Apollo 
11, 14 and 15 sites and on the Lunokhod 1 and 2 rovers. They contribute 10.1%, 10.0%, 
76.9%, 0.2%, and 2.8% of the total observations, respectively. The Apollo 15 retro- 
reflector gives the strongest signal because it is largest. The Lunokhod 1 location was 
only recently determined with sufficient accuracy for ranging after being photographed 
by the Lunar Reconnaissance Orbiter [13]. It gives a stronger return signal than the 
Lunokhod 2 retro-reflector and will be useful in the future. LLR data are archived by 
the ILRS [19]. More information on the ranging stations and technique is in [9]. 

Figure [IK shows the post-fit weighted rms residuals for each year from 1970 to 
2011. The prominent decreases correspond to equipment improvements. Ranges from 
2005-2011 are fit with a 1.8 cm weighted rms residual. The weights depend on the range 
accuracy provided by each station and the analysis software's ability to fit the ranges. 
For the equivalence principle test, LLR analyses are looking for a cosD signature so the 
distribution with D, the mean lunar elongation from the Sun, is important. Figure [T]d 
shows a histogram of the distribution of LLR observations vs. elongation angle D for 
all of the data. The bright Sun prevents observations near new Moon and there are few 
observations near full Moon. Most of the observations near full Moon are less accurate 
early observations. Alone among the modern more accurate stations, the APOLLO 
station provides some ranges near full Moon [12]. The nonuniform distribution with D 
weakens the LLR equivalence principle test. 



4. Data Analysis 

This section presents analysis of the lunar laser ranging data to test the equivalence 
principle. To check consistency, more than one solution is presented. Solutions are 
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made with two different equivalence principle parameters and several different choices 
for other parameters. The solutions are compared and discussed, and residuals after 
fits, the post-fit residuals, are examined for systematics. 

An EP violation can be solved for in two ways. The first is a parameter for 
{mG/mi)E — {rnG/rni)u with a dynamical partial derivative generated from numerical 
integration. The second solves for a coefficient of cos D in the lunar range, a one-term 
representation. The latter approach was used in early papers, but the more sophisticated 
dynamical parameter is used in subsequent publications [28| [23]. Both approaches are 
exercised here to investigate consistency. 

Six equivalence principle solutions are presented in Table [H Each of these cases 
solves for a standard set of Newtonian parameters in addition to one or more equivalence 
principle parameters. Among the Newtonian solution parameters are GmEarth+Moon, 
lunar orbit parameters including semimajor axis. Moon-centered retroreflector 
coordinates, geocentric ranging station coordinates, lunar tidal displacement Love 
number /i2, and parameters associated with the lunar orientation (physical librations) 
and physical properties. The interior model for the physical librations includes a fluid 
core, but the size of that core remains uncertain. Consequently, two different fluid 
core moments were used, 3 x 10~^ and 7 x 10^^ compared to the total lunar moment. 
These cases have designations starting with EP3 and EP7, respectively. The annual 
nutation of the Earth's pole direction in space involves four parameters, the in-phase 
and out-of-phase components of two nutation angles. We normally solve for these four 
nutation parameters in LLR solutions because they are sensitive to the flattening of the 
Earth's core/mantle boundary, and this was done for cases EP31 and EP71. The LLR 
uncertainties for the two in-phase components are about 0.23 milliarcseconds (mas). 
The out-of-phase components have 0.17 mas LLR uncertainties. In these two solutions 
the equivalence principle correlates 0.90 (ignoring signs) with the two in-phase annual 
nutation components, which weakens the EP uncertainty. Consequently, the remaining 
solutions fix the four annual nutation components at the lAU 2000A nutation model 
values given in the 2003 lERS conventions [8J. This is justified because the nutation 
model is compatible with accurate VLBI solutions for terrestrial nutations at the level 
required for the study here. When the annual nutation components are fixed the 
EP uncertainty improves by more than a factor of two as solutions EP32 and EP72 
show. While the first four cases in the table solve for {mG/mj)^ — (^g/''^/)m using the 
numerically integrated partial derivative, the EP73 case solves for coefficients of cos D 
and sin D in range as an alternative. Only the cos D coefficient behaves like the EP. The 
sinD component shows nothing above the noise level. The last case, EP74, provides a 
check of the Eq. (jTj) scaling factor S that relates {mG/mj)^ — ("^g/''^/)m violations to 
the cosD coefficient, as discussed below. The column labeled "converted cosD coef" 
converts {mG/mi)^ — {rnG/mi)M to the coefficient of cosD in radial distance using 
the factor S = —2.94 x 10^'^ mm in Eq. (jTj). The sum of the converted and solution 
cosD coefficients is 0.69 ±3.80 mm, using the 0.9790 correlation, in agreement with the 
estimates of the single EP parameter cases. 
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Table 1. Six solutions for the equivalence principle. 



Solution 


Ann Nut 


imG/mi)E- 


Converted 


cosD coef 


sinD coef 


ID 




{mG/mj)M 


cos D coef 






Units 




X 10-13 


mm 


mm 


mm 


EP31 


solved 


1.03 ±3.39 


-3.03 ±9.97 






EP71 


solved 


0.65 ±3.35 


-1.91 ±9.85 






EP32 


fixed 


0.42 ±1.32 


-1.24 ±3.88 






EP72 


fixed 


0.03 ±1.28 


-0.08 ±3.75 






EP73 


fixed 






0.90 ±3.78 


2.43 ±2.14 


EP74 


fixed 


2.80 ±6.32 


-8.24 ±18.57 


8.93 ±18.50 


2.35 ±2.14 



Solar radiation pressure, like the acceleration from an equivalence principle 
violation, is aligned with the direction away from the Sun and it produces a perturbation 
with the 29.53 d synodic period. This force on the Earth and Moon must be considered 
for the most accurate tests of the equivalence principle. This acceleration is not currently 
modeled in the JPL software. Here we rely on the analysis of Vokrouhlicky [32] who 
considered incident and reflected radiation for both bodies plus thermal radiation from 
the Moon. He finds a solar radiation perturbation of —3.65 ± 0.08 mm cosD in the 
radial coordinate. 

Thermal expansion of the corner cube mounts is another effect worth consideration. 
The peak-to-peak variation of surface temperature at low latitudes on the Moon is 
nearly 300°. The lunar "day" is 29.53 days long. This is the same period as the main 
equivalence principle term so a systematic effect from thermal expansion is expected. 
The phase of the thermal cycle depends on the retro-reflector longitude. If the Apollo 
arrays share the same temperature variations as the surface, then the total vertical 
variation of thermal expansion will be 1 to 2 mm and the cos D-like amplitude would 
be about 2/3 that {2/p for a square wave). The Apollo arrays make up 97% of the 
data so the range effect could be as much as 1.0 mm cosD, but will be less if the array 
temperature variations are less than the lunar surface variations. 

For the sum of the thermal and solar radiation pressure effect we use —3.0 ±0.5 mm 
cosD. The correction to Table [1] has opposite sign; this correction is used to construct 
Table [2] For the EP74 case the partitioning of the correction for the two highly correlated 
EP parameters is not clear and that case is not given in Table O However, the 
sum of those two EP parameters would give 3.69 ± 3.83 mm cosD, equivalent to 
(-1.26 ± 1.30) X 10-13 for (mG/m/)E - (mG/m/)M. 

The EP solution parameters in Tables [H and [2] are within their uncertainties for 
all cases except the EP73 case in Table |2l and that value is just slightly larger. The 
EP73 coefficient of cos D and the sum of terms for EP74 are about 1 mm larger than 
the value from the converted {rnG/rni)E — {jnGlf^i)M parameter of the EP72 case. The 
uncertainties are compatible. Solutions with the two lunar interior models differ by 
about 1 mm with the smaller lunar core giving smaller cosD coefficients. Solutions 
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Table 2. Five solutions for the equivalence principle corrected for solar radiation 
pressure and thermal expansion. 



Solution 
ID 


Ann Nut 


{mG/mi)M 


Converted 
cosD coef 


cosD coef 


sinD coef 


Units 




X 10^^3 


mm 


mm 


mm 



EP31 solved 

EP71 solved 

EP32 fixed 

EP72 fixed 

EP73 fixed 



0.01 ±3.40 
-0.37 ±3.35 
-0.60 ±1.33 
-0.99 ±1.29 



-0.03 ±9.98 
1.09 ±9.86 
1.76 ±3.91 
2.92 ±3.78 



3.90 ±3.81 2.43 ±2.14 



with different EP parameters are compatible. There is no evidence for a violation of the 
EP from any of the solutions. 

The EP32 and EP72 solutions are the preferred cases since they use numerically 
integrated partial derivatives for {mG/'mj)^ — {fnG/mi)u- The average of those two 
corrected solutions is {mG/mi)^ — {jnG/mi)y[ = (—0.8 ±1.3) x 10"^^. Some larger 
correlations from the EP72 solution follow. The correlation of {mG/mi)^ — {fnG/rni)M 
with both GmEarth+Moon and osculating semimajor axis, at the 1969 epoch of the 
integration, is 0.31. Gm and mean semimajor axis are tightly connected through 
Kepler's third law since the mean motion is very well determined. The product of 
mean semimajor axis and mean eccentricity is well determined and the correlation of 
{mG/mi)-E — {jnG/rni)M with osculating eccentricity is 0.29. The displacement Love 
number /i2 is correlated 0.28 with {mG/mj)-^ — {mG/mj)y[- AH of these correlations are 
modest. 

The difference in uncertainty between the sin D and cos D components of both the 
EP73 and EP74 solutions is due to the nonuniform distribution of observations with 
respect to D, as illustrated in Fig. [lb. The sinD coefficient is well determined from 
observations near first and last quarter Moon where sinD has its extreme values of +1 
and —1, respectively. The cosD coefficient is weakened by the absence of data close to 
new Moon and its scarcity near full Moon, where cos D has its extreme values. 

The EP74 case, solving for {mG/mi)-^ — {mG/mi)^ along with cosD and sinD 
coefficients, is instructive. The correlation between the {mG/rni)E — {rnG/mi)u and 
cosD parameters is 0.9790 so the two quantities are nearly equivalent, as expected. 
The uncertainty for each of the two cos D coefficients increases by a factor of five in 
the joint solution, but the sum is 0.69 ± 3.80 mm, which corrects to 3.69 ± 3.83 mm, in 
agreement with the other solutions. The solution is not singular, so the solution has some 
ability to distinguish between the two formulations. The integrated partial derivative 
implicitly includes terms at frequencies other than the D argument and it will also have 
some sensitivity to the EP influence on lunar orbital longitude. The EP perturbation on 
lunar orbital longitude is about twice the size of the radial component and it depends on 
sinD. The ratio of Earth radius to lunar semimajor axis is R^/a ~ 1/60.3, the parallax 
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Table 3. Theoretical and numerical results for factors 5* and Cp. 

Source Type S Cq 

Nordtvedt (1995) [16j analytical -2.9 x lO^^ ^^ 12.90 m 
Damour & Vokrouhlicky (1996) [5] analytical -2.9427 x lO^^ mm 13.10 m 
Nordtvedt & Vokrouhlicky (1997) [TS] analytical -2.943 x lO^^ ^^ x3_;l0 m 

This paper numerical —2.992 x 10^^ mm 13.31 m 

is ~ 1°, so the sinD longitude component projects into range at the few percent level. 

In the radial cosD expressions the coefficients of the EP equations ([5]) and (jTj) 
have an S factor, while the strong EP Cq factor of Eq. ([8]) is the product of S and the 
difference in Earth and Moon self-energies. Three theoretical values for S are given in 
Table [3] while the Co values are based on the self-energy difference given in Eq. ([6]). The 
uncertainties in the EP74 solution can be used to check the theoretical computation of 
the coefficients S and Cq. Given the high correlation between the {jnG / mi)-E, — {rnG / mi)^ 
and cosD parameters, a ffist approximation oi S = —2.93 x 10^^ mm is given by the 
ratio of uncertainties, and our knowledge that 5* must be negative. A more sophisticated 
estimate oi S = —2.992 x 10^'^ mm comes from computing the slope of the axis of the 
uncertainty ellipse for the two parameters Co and {mc / rni)-^. — {jnG / mi)yi. Using Eq. ([6]) 
for the difference in self energies of the Earth and Moon, the two preceding S values 
give Ar = 13.0 m rj cos D and Ar = 13.3 m rjcosD, respectively. As explained in the 
preceding paragraph, the LLR solution for {mG/mj)-^ — {mG/rnj)M involves more than a 
cos D radial expression so we expect a few percent uncertainty in the numerical values. 

Our earlier paper [30] discusses the importance to the LLR equivalence principle 
solutions of parameters that affect the mean distance of the lunar retroreflectors 
including GmEarth+Moon and lunar displacement Love number /i2. Both were solution 
parameters in all cases of Tables [1] and |2l 

Six solutions presented in this analysis section have tested the equivalence principle. 
They do not show evidence for a significant violation of the equivalence principle. 

5. Future Prospects 

In the fits, the data is weighted using both uncertainties assigned by the ranging stations 
and a contribution from the analysis software. The latter depends on the post-fit rms 
residuals by station and time. Modern high-quality ranges have smaller instrumental 
uncertainties than the analysis software's ability to fit over decades-long time spans. 
There have been software improvements since [281 EH] but further improvements are 
needed. Most obvious in the highest accuracy residuals are systematic signatures that 
are different for different retro-reflectors. This indicates that there are unmodeled 
variations in the physical librations. The structure and geophysical properties of the 
deep lunar interior are imperfectly known. Further software improvements will advance 
the LLR equivalence principle tests. These upgrades are planned. 
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The accuracy of the annual nutation will become more of a concern for future 
improved LLR tests of the equivalence principle. Four of the solutions given in this 
paper fix the values of the four annual nutation coefficients and the EP uncertainties do 
not account for uncertainty in the annual nutation. It is estimated that 0.05 milliarcsec 
uncertainties in the two in-phase annual nutation coefficients would increase the LLR 
EP uncertainties by about 10%. 

Future ranging devices on the Moon might take two forms, namely passive retro- 
reflectors and active transponders. The advantages of passive retro-refiector arrays are 
their long life and simplicity. The disadvantages are the weak returned signal and the 
spread of the reflected pulse arising from lunar orientation with respect to the direction 
toward the Earth, which can alter the incoming direction up to 10° with respect to 
the retro-reflector face. A single large corner cube would not have this problem |2]. 
Additional ranging devices on the Moon would have beneflts for fundamental physics, 
lunar science, control networks for surface mapping, and navigation [231 I2H [25] . 

6. Summary 

If the ratio of gravitational mass to inertial mass is not invariant, then there would be 
profound consequences for gravitation. Such a violation of the EP would affect how 
bodies move under the influence of gravity. The EP is not violated in general theory 
of relativity, but violations are expected for many alternative theories of gravitation. 
Consequently, tests of the EP are important to the search for a new theory of gravity. 

The equivalence principle (EP) is considered in its two forms; the weak equivalence 
principle (WEP) is sensitive to composition while the strong equivalence principle (SEP) 
considers possible sensitivity to the gravitational energy of a body. The main sensitivity 
of the lunar orbit to the equivalence principle comes from the acceleration of the 
Earth and Moon by the Sun. Any difference in those accelerations due to a failure 
of the equivalence principle causes an anomalous term in the lunar range with the 
29.53 d synodic period. The amplitude would be proportional to the difference in the 
gravitational to inertial mass ratios for Earth and Moon. Thus, lunar laser ranging is 
sensitive to a failure of the equivalence principle due to either the WEP or the SEP. In 
the case of the SEP, any violation of the equivalence principle can be related to a linear 
combination of the parameterized post-Newtonian parameters /3 and 7. 

Section [3] discusses the data and its distribution with respect to the new-full-new 
Moon cycle. The evolution of the data from decimeter to centimeter quality flts is 
illustrated. For the LLR equivalence principle tests, selection with phase of the Moon 
is an important consideration. 

An accurate model and analysis effort is needed to exploit the lunar laser range 
data to its full capability. The model is the basis for the computer code that processes 
the range data. Further modeling efforts will be necessary to fully beneflt from range 
data of millimeter quality. Two small effects for future modeling are thermal expansion 
and solar radiation pressure. An improved model of the deep lunar interior is needed 
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for lunar orientation. 

Solutions testing the EP are given in Sec. IH Several approaches to the solutions 
are used as checks. The equivalence principle solution parameter can be either a ratio 
of gravitational to inertial masses or as a coefficient of a synodic term in the range 
equation. The results are compatible in value and uncertainty. In all, six equivalence 
principle solutions are presented in Table [1] and five are corrected and carried forward 
into Table |2j The analysis of the LLR data does not show significant evidence for a 
violation of the equivalence principle compared to its uncertainty. The final result for 
(mG/m,)E - {mG/mj)u is (-0.8 ± 1.3) x lO-^^. 

The accuracy of the LLR test of the EP can be improved in the future. Adding 
new LLR data will help. The analysis software needs to be improved so that the most 
accurate data, including the very accurate APOLLO data, can be more fully exploited. 
Future landers on the Moon should carry new retro-refiectors with improved design. 

The lunar laser ranging results for the EP are consistent with the assumptions 
of Einstein's general theory of relativity. It is remarkable that general relativity has 
survived a century of testing. Each new significant improvement in accuracy is unknown 
territory and that is reason for future tests of the equivalence principle. 
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